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The  aim  of  this  paper  is  to  develop  and  manufacture  a  thermoelectric  structure  operating 
in  Peltier  mode  with  non-standard  materials,  using  a  Finite  Element  Model.  Thus,  both  the 
reliability  of  the  Finite  Element  Model  and  the  correct  development  and  design  of  the  new 
thermoelectric  structure  are  ascertained,  achieving  a  dual  goal  of  our  research  group.  This 
work  provides  a  very  good  correlation  between  simulated  and  experimental  results,  and 
corroborates  that  the  material’s  performance  has  a  decisive  influence  on  the  thermo¬ 
electric  properties  of  the  structure  in  the  temperature  range  of  study  considered. 
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Nomenclature 

Rua>b  Newton— Raphson  electrical  residual  vector 

RTA)B  Newton— Raphson  thermal  residual  vector 

Kw.vt.tv.vt  Derivatives  of  the  residual  vectors 
dVA>B  Unknown  regarding  to  voltage  in  the  Finite 
Element  equations  systems 

dTA  B  Unknown  regarding  to  temperature  in  the  Finite 

Element  equations  systems 
Na> b  Shape  spatial  functions 

jh  Discretized  electrical  current  flux 

Vh  Discretized  voltage 

Th  Discretized  temperature 

Ipeiiet  Electric  current  through  a  pellet  (A) 

IT  Electric  current  through  a  thermoelectric  module 

(A) 

Mn  n-material 

Mpi  p-material  type  1 

Mpi  p-material  type  2 

Qt  Total  heat  flux  of  a  thermoelectric  module  (W) 
Qpeiiet  Pellet  heat  flux  (W) 

Tc  Cold  side  temperature  (K) 

Th  Hot  side  temperature  (K) 

Vpeiiet  Voltage  of  a  pellet  (V) 

VT  Voltage  of  a  thermoelectric  module  (V) 

ZT  Dimensionless  Figure  of  Merit 

Greek  symbols 

a  Seebeck  coefficient  (jiV  K-1) 

k  Thermal  conductivity  (mW  cm-1  K-1) 


a  Electrical  conductivity  (Q-1  cm-1) 

r  Closed  integration  surface 

Q  Contour  volume 

rjv  Set  of  nodes  belong  to  rv 

r]T  Set  of  nodes  belong  to  rT 

qh  Discretized  thermal  flux 

ah  Discretized  Seebeck  coefficient 

Kh  Discretized  thermal  conductivity 

Discretized  electrical  conductivity 

Formulas 

Bi2Te3  Bismuth  telluride 
Sbl3  Antimony  triiodide 
Sb2Se3  Antimony  selenide 
Sb2Te3  Antimony  telluride 

Subscripts 

A,  B  Nodes  in  Finite  Element  formulation 
T  Regarding  to  temperature  or  heat  transfer 

V  Regarding  to  electrical  phenomena 

i  Iteration  number  in  the  Newton— Raphson 

method 

C  Cold  side 

Ce  Ceramics 

Cu  Copper 

H  Hot  side 

S  Welding 

Superscript 

h  Discretized  variable 


1.  Introduction 

Since  the  discovery  of  thermoelectric  phenomena  two  cen¬ 
turies  ago  (Seebeck  effect,  Peltier  effect,  and  Thomson  effect), 
their  application  in  both  heat  pump  (Peltier  effect)  and 
electrical  power  generation  (Seebeck  effect)  has  experienced 
an  important  increase  during  the  last  60  years. 

Thermoelectric  technology  is  currently  spreading  across 
many  application  areas  such  as  power  generation  from 
waste  heat  (Kajikawa,  2006),  thermoelectric  systems  in  the 
automotive  industry  (Matsubara  and  Matsuura,  2006),  as  well 
as  being  applied  to  space  missions  (Abelson,  2006)  and  for 
cooling  electronic  components  (Semenyuk,  2006). 

Therefore,  this  article  describes  the  development  of  a 
thermoelectric  module  made  of  non-standard  p-n  semi¬ 
conductors.  These  materials  are  compounds  of  type 
(Bi2Te3)a_x_y  (Sb2Te3)x  (Sb2Se3),  and  show  good  thermoelectric 
properties  as  n-type  and  p-type  materials  within  a  tempera¬ 
ture  range  from  270  K  to  450  K.  These  p-n  semiconductors 
have  been  obtained  through  the  Bridgman-Stockbarger  tech¬ 
nique  (Kuznetsov  et  al.,  2002).  It  should  be  noted  that  the 
module  suggested  in  this  paper  is  intended  to  perform  within 
a  food  cooling  system  and  a  cigar  humidor,  due  to  its  low 
power  consumption  and  its  high  refrigerating  power. 

Nowadays,  authors  such  as  Cai  et  al.  (2011)  have  used 
analytical  models  to  develop  thermoelectric  structures, 
obtaining  optimal  results.  Other  authors  have  used  numerical 


methods  to  simulate  the  behaviour  of  thermoelectric  modules 
(Yu  and  Zhao,  2007).  In  this  paper  a  discrete  model  based  on 
Finite  Elements  was  used.  This  model  has  been  already  tested 
in  commercial  structures  and  other  internal  developments 
through  specific  software  (Perez-Aparicio  et  al.,  2007,  2012). 
The  results  have  been  satisfactory,  thus  validating  its  utility. 


2.  Thermoelectric  properties  of  materials 
considered 

We  have  developed,  characterized  and  tested  thermoelectric 
plates  with  pellets  made  of  materials  generated  and  charac¬ 
terized  within  the  research  group.  The  materials  used  are  the 
following: 

-  MP1:  p-type  material  (Bi2Te3)0.3  (Sb2Te3)0.67  (Sb2Se3)0.o3  doped 
with  (Sbl3,  Te)  type  1. 

-  MP2:  p-type  material  (Bi2Te3)0.3  (Sb2Te3)0.67  (Sb2Se3)0.o3  doped 
with  (Sbl3,  Te)  type  2. 

-  Mn:  n-type  material  (Bi2Te3)0.s  (Sb2Te3)0.i  (Sb2Se3)0.i  doped 
with  (Sbl3). 

These  materials  can  form  cells  of  two  types:  one  with  MP1 
and  Mn  pairs,  and  another  one  with  MP2  and  MN  pairs.  This 
paper  shows  the  development  of  a  structure  with  MP1  and  MN 
materials.  The  thermoelectric  properties  at  room  temperature 
(300  K)  of  these  materials  are  shown  in  Table  1. 


1572 


INTERNATIONAL  JOURNAL  OF  REFRIGERATION  36  (2013)  I57O-I575 


Table  1  -  Thermoelectric  properties  of  materials  used  (T  =  300  K). 


Material 

Type 

Seebeck  coefficient 

Electrical  conductivity 

Thermal  conductivity 

Figure  of  Merit  ZT 

MP1 

P 

160 

1900 

19.3 

2.4 

Mp2 

P 

182 

1290 

15.2 

2.8 

Mn 

n 

-176 

1710 

19 

2.8 

The  MP1  sample  has  a  higher  carrier  concentration  compared  to  the  MP2  sample.  Therefore,  MP1  has  higher  electrical  and  thermal  conductivity. 
Both  samples  show  a  thermal  conductivity  lower  than  Bi2Te3  (Ivanova  et  al.,  2002). 


3.  Design  and  simulation  model 


dv  =  {dvB}  Be  77  —  ?7V 


(9) 


The  Finite  Element  Model  used  allowed  us  to  find  a  spatial 
distribution  of  electric  potential  and  temperature  in  a  volume 
Q  contained  in  a  closed  surface  r,  with  specific  electrical  and 
thermal  boundary  conditions  (Virjoghe  et  al.,  2010).  The  model 
is  based  on  second-order  differential  equations,  which  are 
simplified  (first-order  integral)  and  discretized  for  application 
in  the  Finite  Element  theory,  thus  obtaining  a  system  of  non¬ 
linear  equations  which  we  solved  with  the  Newton— Raphson 
method  (Gavela  and  Perez-Aparicio,  1998),  as  shown  below: 


Kwab  —  - 


9RVa 

0dyB 


Kvtab 


0RVa 

ddTB 


Q 


A,  Be  77  -  77V 


Ae  77  —  77v,  Be  77  —  77t 


(1) 

(2) 


Ktvab 


9Rta 

9dyB 


J  NAjhVNBdQ 

Q 


Ae  77 


—  77t,  Ber]  —  riv 


(3) 


Q  Q 

(4) 

From  Equation  (4),  the  value  of  the  derivatives  can  be  ob¬ 
tained  according  to  the  intrinsic  characteristics  of  the 
semiconductor: 


a f 

ddVB 


=  -<rhVN  j 


(5) 


a f 

ddjB 


=  -cth(ThVNB  -  -^r<rhVThNB 
dT 


^(VVh  +  ahVTh)NB 


(6) 


aql  =  aHT^ 

ddyB  ddyB 


(7) 


dT  =  {dTB}  Be  77  —  77t  (10) 

Then  we  can  apply  the  Newton— Raphson  iterative  process 
in  order  to  find  the  values  for  dv  and  dT.  The  process  can  be 
described  in  matrix  form  as  follows,  where  i  is  the  corre¬ 
sponding  iteration: 


In  Equations  (5)— (8),  this  three-dimensional  model  can 
replace  the  intrinsic  characteristics  of  the  MP1  and  MN  semi¬ 
conductors,  used  in  the  formation  of  p-n  pairs.  Properties  of 
the  remaining  elements  in  the  thermoelectric  system,  such  as 
copper,  aluminium  oxide,  and  welding,  are  considered. 

Copper  shows  values  of  aCu  =  1-84  pV  K-1, 

o-Cu  =  5.81- 107  Q-1  m-1  and  kCu  =  386  W  m_1,  with  electrical  and 
thermal  conductivities  assumed  to  be  temperature- 
independent.  In  the  case  of  aluminium  oxide,  with  which  ce¬ 
ramics  is  made,  we  obtain  a  value  of  KCe  =  35.3  W  m-1,  also 
assumed  to  be  temperature-independent.  Finally,  in  the  joint 
welding  (SnPb)  between  the  thermoelements,  we  obtain 
o s  =  17  -104  Q-1  m_1  and  ks  =  48  W  m-1,  also  with  constant 
properties.  The  values  of  the  Seebeck  coefficients  of  p-n  mate¬ 
rials  and  the  remaining  elements  of  the  thermoelectric  system 
have  been  obtained  experimentally  by  performing  a  Seebeck 
Scanning  Microprobe  (Fig.  1),  as  detailed  in  Platzek  et  al.  (2005). 

In  the  simulation,  it  is  assumed  that  all  pairs  of  thermo¬ 
couples  behave  the  same  way.  This  depends  mainly  on  the 
temperature  distribution  between  the  two  faces  of  the  pellet 
(Antonova  and  Looman,  2005).  If  the  temperature  on  the  faces 
is  homogeneous,  all  pairs  of  thermocouples  work  under  the 
same  conditions.  Therefore,  including  only  one  pair  of  ther¬ 
mocouples  is  enough,  considering  the  top  as  the  cold  face,  and 
the  bottom  as  the  hot  face.  Electric  current  enters  the  model 
through  the  bottom  left  side  of  the  copper  conductor,  and 
leaves  it  through  the  bottom  right  side. 

It  was  decided  to  build  a  board  of  49  pairs  of  thermocou¬ 
ples,  with  pellets  made  of  the  materials  considered,  MP1  and 
Mn.  These  pellets  have  a  size  of  3  mm  x  3  mm  x  3  mm.  The 
total  sizing  of  the  thermoelectric  system  with  distribution  of 
the  49  pairs  (98  pellets)  is  shown  in  Figs.  2  and  3. 

3.1.  Boundary  conditions 


|f-  =  +  «hjhNB  -  KhVNB  +  ^Thjh Nb  -  ^pVThNB 

ddT]}  ddrB 


dT  J  D  dT 
If  dv  and  dT  are  vectors  containing  the  unknowns: 


(8) 


The  boundary  conditions  have  been  established  as  follows: 

-It  is  assumed  that  the  entire  contour  is  open.  Therefore,  the 
contour  is  not  thermally  insulated  and  the  ceramic  faces 


INTE  RNATIONAL  JOURNAL  OF  REFRIGERATION  36  (2013)  I57O-I575 


1573 


Fig.  1  —  Scanning  Seebeck  microprobe  carried  out  on  a  MP1 
pellet. 


(upper  and  lower)  are  considered  at  homogeneous  tempera¬ 
tures  Th  and  Tc.  Moreover,  there  is  convection  in  the  ther¬ 
mocouples.  After  the  transient  period  Th  and  Tc  tend  to  be 
constant. 

-It  is  assumed  that  the  entire  contour  is  electrically  insulated, 
except  for  the  ends  of  the  copper  conductors  through  which 
electric  current  enters  and  leaves. 

The  model  supports  three  input  parameters:  the  electric 
current  through  the  module,  the  temperature  of  the  hot  face, 


and  the  temperature  of  the  cold  face.  The  model  calculates 
the  voltage  and  temperature  in  domain  points,  obtaining 
both  current  and  heat  flux.  From  this  spatial  distribution,  we 
can  also  obtain  important  operating  parameters  of  the  pel¬ 
lets,  such  as  the  voltage  applied  to  the  couple,  the  heat 
extracted  from  the  cold  face  and  the  heat  dissipated  on  the 
hot  face. 

This  last  parameter  is  not  independent  of  the  rest,  since 
the  heat  dissipated  on  the  hot  face  must  equal  the  heat 
absorbed  by  the  cold  face  plus  the  electrical  power 
consumed.  However,  the  heat  dissipated  on  the  hot  face 
obtained  by  the  model  confirms  that  everything  works 
properly.  In  order  to  relate  the  operating  conditions  of  the 
pellet  to  the  parameters  of  the  model,  we  will  take  the  con¬ 
siderations  described  below. 

3.2.  Input  parameters 

Electric  current.  All  thermoelectric  couples  are  connected  in 
series  and  supplied  by  constant  voltage,  which  means  that  the 
total  electrical  power  consumed  by  the  cell  flows  through  the 
couples.  This  model  includes  49  couples,  so  the  electric  cur¬ 
rent  applied  to  the  model  will  be  the  same  that  flows  through  a 
pellet,  that  is  It  =  Icoupie  =  Ipeiiet- 

3.3.  Output  parameters 

Voltage.  The  total  voltage  applied  to  the  thermoelectric  sys¬ 
tem  is  the  sum  of  all  couples,  that  is,  VT  =  98  x  Vpellet.  To 
maintain  the  value  of  the  electric  current  Ipeiiet  it  is  necessary 
to  supply  the  system  with  a  voltage  VT. 

Heat  flux.  The  heat  exchanged  between  the  plate  and  the 
faces  is  the  amount  of  heat  for  every  couple,  that  is, 
Qt  =  98  x  Qpgiieti  assuming  no  losses  in  all  associated  elements 
(ceramics,  welding,  etc.). 

The  results  obtained  in  the  simulations,  taking  into  ac¬ 
count  the  complete  structure  parameters,  are  described  in 
detail  in  the  following  sections. 


■  p-type  semiconductor 

■  n-type  semiconductor 

■  solder 
ceramic  plate 


Fig.  2  -  Design  scheme  of  the  thermoelectric  module. 
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Fig.  3  —  Thermal  structure  developed. 

4.  Structure  study  operating  in  Peltier  mode 

4.1.  Results  obtained  in  the  simulation 

Considering  that  the  thermoelectric  structure  is  working  in 
Peltier  mode,  the  following  results  were  obtained  with  250 
differential  elements  in  the  iteration.  In  these  simulations  we 


Fig.  5  -  Experimental  measurements  of  the  temperature 
evolution  on  the  hot  and  cold  faces  for  different  current 
values. 


can  see  the  theoretical  behaviour  of  the  plate,  which  can  then 
be  compared  to  the  measurements  carried  out  on  the  devel¬ 
oped  structure.  The  theoretical  development  of  temperature 
in  the  time  domain  for  different  values  of  input  current  (1A, 
1.5A,  2A,  2.5A  and  3A)  for  100  s  is  shown  in  Fig.  4. 

4.2.  Experimental  results 

The  measurements  on  the  designed  and  developed  structure 
were  performed  with  the  same  values  of  electrical  current 


Fig.  4  —  Simulation  of  the  temperature  evolution  on  the  hot  Fig.  6  -  Difference  between  experimental  and  theoretical 
and  cold  faces  for  different  current  values.  values  on  the  cold  face  for  different  current  values. 
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Time  [s] 


Fig.  7  -  Difference  between  experimental  and  theoretical 
values  on  the  hot  face  for  different  current  values. 

(1A,  1.5A,  2 A,  2.5A,  and  3A).  Temperatures  TH  and  Tc  were 
obtained  for  each  current  value,  as  shown  in  Fig.  5. 

As  shown  in  Figs.  6  and  7,  although  the  simulation 
response  is  slightly  faster  while  booting  the  system  than  in 
the  experimental  module,  the  difference  between  the 
experimental  and  theoretical  results  decreases  over  time. 
Therefore,  the  difference  is  minimal  and  the  model  behaves 
correctly. 


5.  Conclusions 

In  this  paper  we  have  characterized  a  thermoelectric  structure 
with  p-type  thermocouples  formed  by  (Bi2Te3)0.3  (Sb2Te3)0.67 
(Sb2Se3)0.o3  doped  with  (SbI3Te)  type  1,  in  the  temperature 
range  from  270  K  to  450  K.  After  the  different  simulations  and 
tests  have  been  executed,  we  can  deduce  a  direct  relationship 
between  the  thermoelectric  properties  of  the  module  and  the 
properties  of  the  constituent  material  of  the  pellets.  There¬ 
fore,  we  can  confirm  the  expectations  of  our  previous  studies, 
which  allowed  us  to  sense  the  influence  of  the  pellet  proper¬ 
ties  on  the  behaviour  of  the  thermoelectric  module. 

Comparisons  between  simulation  results  and  experi¬ 
mental  results  show  a  very  good  correlation,  especially 
regarding  the  heat  extracted  from  the  cold  face.  Furthermore, 
it  should  be  noted  that  small  deviations  in  temperature  are 
largely  due  to  the  not  well-known  values  of  the  actual  con¬ 
ductivity  and  thickness  of  the  copper  junctions  of  the  ther¬ 
mocouple  implemented  during  the  development  of  the 
module.  Moreover,  it  should  be  pointed  out  that  deviations 
obtained  by  Cai  et  al.  (2011)  and  Yu  and  Zhao  (2007)  with  their 


modules  are  in  the  same  order  of  magnitude  as  those  obtained 
with  the  model  proposed  in  this  paper. 

Thus,  in  addition  to  characterising  the  thermoelectric 
structure,  the  reliability  of  the  Finite  Element  Model  used  has 
been  verified,  since  the  theoretical  values  obtained  are  very 
similar  to  the  experimental  results.  Furthermore,  as  discussed 
above,  the  influence  of  thermoelectric  materials  on  the  mod¬ 
ule  behaviour  in  the  considered  temperature  range  of  study 
has  been  confirmed  to  be  decisive. 
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